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Exact diagonalization study of a rotating dilute Bose-Einstein condensate reveals that as the first 
vortex enters the system the degeneracy of the low-energy yrast spectrum is lifted and a large energy 
gap emerges. As more vortices enter with faster rotation, the energy gap decreases towards zero, 
but eventually the spectrum exhibits a rotonlike structure above the v = 1/2 Laughlin state without 
having a phonon branch despite the short-range nature of the interaction. 
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Rotating dilute Bose-Einstein condensates (BECs) 
have attracted considerable interest in recent years [ij, 
EJ When the single-particle energy- level spacing 
is much larger than the interaction energy per particle, 
the system becomes highly degenerate with increasing 
the angular momentum (AM) L. Interactions then play 
a pivotal role in lifting the degeneracy and determining 
the many-body ground state |5(, in close analogy with 
the physics of the fractional quantum Hall effect 0, • 
Mottelson on the other hand, pointed out the rel- 
evance of the yrast line, which traces the lowest-lying 
states of the system as a function of the AM and has 
served as a key concept in nuclear physics 0, to a ro- 
tating dilute BEC. In view of recent development of Fes- 
hbach techniques and experimental achievements of 
fast rotating BECs [l], [2J, |3( , such an extreme but highly 
interdisciplinary arena seems to be within experimental 
reach. 

Bertsch and Papcnbrock numerically found that the 
yrast line depends linearly on L for L < N, where N 
is the number of bosons and showed that high- 

lying states above the yrast line are dominated by single- 
particle excitations |12|. Of particular interest is the 
yrast spectrum for low-lying states because it gives us in- 
sights about how the many-body wave function responds 
as vortices enter the system. Mottelson pointed out that 
for L -C N low-lying states are quasi-degenerate and 
dominated by collective multipolar excitations ||. An 
exact diagonalization study [13| and subsequent analyti- 
cal studies [Til Hil fl6| have demonstrated the existence 
of the quasi-degenerate yrast spectrum that is dominated 
by interactions between octupole modes. The purpose of 
this Letter is to report the results of our exact diagonal- 
ization study on faster rotating BECs. Our primary find- 
ings are two distinct kinds of energy gaps: one above the 
single-vortex state and the other associated with vortex- 
antivortex pair excitations above the v = 1/2 Laughlin 
state. We also find that higher excited states with L < N 
form linear energy bands. 

Consider a weakly interacting A^-boson system trapped 



in an axisymmetric parabolic confining potential. The 
potential is assumed to be isotropic in the radial di- 
rection but in the axial direction it is assumed to be 
so tightly confined that all particles occupy the lowest- 
energy state in that direction. The problem is thus 
essentially two-dimensional with projected AM, L, on 
the symmetry axis being conserved. The many-body 
Hamiltonian of this system can be written as H = 
J2i hi + V, where hi = -h 2 Vf/2M + Mu 2 r 2 /2 denotes 
the single-particle Hamiltonian for the z-th particle and 
V = (Airhg/Muj) J2i<j^( r i ~ r j) describes the contact 
interaction between bosons having mass M. Here g 
characterizes the strength of interaction between bosons, 
g/Tiuj — a s /(27r?i/Mu) z ) 1 ' 2 , uj and uj z are trap frequen- 
cies in the radial and axial directions, respectively, and 
a s is the s-wave scattering length. Throughout this Let- 
ter we shall assume weak repulsive interactions such that 
huu > gN > 0. 

The single-particle state is characterized by the ra- 
dial quantum number n and the AM quantum num- 
ber to, with the corresponding eigenenergy given by 
hu>(2n + \m\ + 1). For a given AM L > 0, it is en- 
ergetically favorable to assume n — and let all par- 
ticles have nonnegative AM, i.e., to > 0. In this 
lowest-Landau-level approximation, the noninteracting 
part of the Hamiltonian contributes a constant term 
Uu)(L + N), and we shall henceforth ignore this constant 
part and focus on the interaction Hamiltonian. It is con- 
venient to take as a basis set of states single-particle 
states described by <j) m (z) = (z m / 'V Vm!) exp(— \z\ 2 /2), 
where z = x + iy and the lengths are measured in 
units of (Ti/Moj) 1 / 2 . Expanding the field operator as 
^( z ) — Em'™^ra( z )' t ne second-quantized form of 
the contact interaction Hamiltonian is written as V = 
9Y, mi ,..., mi V rnirn 2m3mi bl ni bl l2 b m3 b m4 , where the ma- 
trix elements V mim2 m 3mi are given by V mi m 2 m 3m4 , = 
<W+m 2 ,m 3 + m4 + m 2 )l/(2 mi+m2 V"ii!TO 2 !m 3 !TO4!). 
Our task thus reduces to finding the energy spectrum 
of V under the restrictions of a fixed particle number 
E m b^bm = N and a fixed AM £) m mb^bm = L. 
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Figure ^ shows the yrast spectrum for < L < 35 
with N = 25, where the energy is measured in units of 
gN from the yrast line E = gN(N -l-L/2) [Hj]. For 
L <C N, there are quasi-degenerate excited states aris- 
ing from pairwise repulsive interaction between octupole 
modes yUj. However, the excitation energies are very 
small, of the order of g. This quasi-degeneracy may be 
regarded as a precursor for spontaneous symmetry break- 
ing of the axisymmetry associated with the entrance of 
the first vortex. Remarkably, the quasi-degeneracy for 
L <C N is lifted with increasing L and that a large en- 
ergy gap of the order of gN appears as the first vortex 
enters the system. The energy gap implies that the many- 
body wave function responds to external rotation in such 
a manner that a single vortex state becomes stabilized. 

The emergence of the energy gap can be inferred 
from the fact that excitations associated with the single- 
particle state 4> rn (m > 3) become increasingly massive, 
as the AM per particle L /N approaches unity. As shown 
in Table J] the excitation energy increases by an amount 
of the order of gN when the condensate develops from 
the nonvortex state at L = to the single-vortex state 
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FIG. 1: The yrast spectrum for < L < 35 with N = 25. As 
the angular momentum per particle L/N approaches unity, 
excitations associated with single-particle states (f> m (m > 3) 
become increasingly more massive. The amount of increase 
in the excitation energy is given in Table [I] Among rotational 
bands stabilized by these energy gaps, three bands are shown 
by linking the collective excitations, Q\\L = N) (solid curve), 
Q\\L = N - 1) (dashed curve), and Q\\L = N — 2) (dotted 
curve), for A < 7. Those states that include excitations of 
the center-of-mass motion are omitted except for the states, 
Qi\L = 0) and Q\\L = N) (shown by squares). The energy is 
measured from the yrast line E = gN(N — 1 — L/2) in units 
of gN. 



at L = N, and an increase in excitation energy is larger 
for smaller m. In particular, excitations of the m = 3 
mode, which are nearly gapless for L/N -C 1, become in- 
creasingly massive as the vortex enters the system. The 
energy gap continues to increase till L = N + m — 1 
for each m. We also note that the energy bands (shown 
in Fig. ^ as dash-dotted lines) for m = 4, 5, 6, • • • are 
almost linear. When L/N > 1, a regular structure of ro- 
tational bands is stabilized by the energy gaps associated 
with the single-particle states <p m . In Figure ^ we show 
three rotational bands by linking the collective excita- 
tions Q\\L = N), Q X \L = N - 1), and Q\\L = N - 2) 
(A < 7) with solid, dashed, and dotted curves, respec- 
tively. Here Q\ = (1/ViVA!) J2p=i z p describes the col- 
lective multipolar excitation of order A @, and \L = N), 
\L = N — 1), and \L = N — 2) denotes the yrast states 
for L = N, L = N - 1, and L = N — 2, respectively. 

The regularity of rotational bands for L/N > 1 is more 
clearly seen in Fig.0 where the energy is measured from 
the line E = E L=N - 5gN(L- N)/32 As in Fig.ffl 
the collective excitations, Q\\L = N), Q\\L = N — 1), 
and Q\\L = N — 2) (A < 7), are linked by solid, dashed, 
and dotted curves, respectively. Those states that con- 
tain excitations of the center-of-mass motion are not 
shown except for the states Q\\L = N), Q\\L = N — 1), 
and Qi\L = N — 2) (shown by squares). Other low-lying 
excitations can be understood as (multiple) excitations 
of these collective modes. 

When L > N(N — 1), the system can respond in 
such a manner that the interaction energy becomes zero. 
In fact, the v = 1/2 Laughlin state |l8|, ^ r i/=i/2 = 

Ilpe- |Zp|2/2 Yli<j(^ - z j) 2 , appears at L = N(N - 1) 
as the first zero-interaction-energy state; the contact in- 
teraction does not affect this state because the minimum 
of the relative AM between particles is two for this state. 
In such a high AM regime, we find a remarkable energy 
gap to appear above the v — 1/2 Laughlin state as shown 
in Fig. □ 

To investigate long-wavelength collective excitations 
while keeping the average density unaltered (i.e., neutral 
excitations), it is convenient to use the spherical geome- 
try |l9| rather than the disk geometry which is used thus 
far. The two geometries are related with each other by 
the stereographic mapping. The z-component of the total 
AM in the spherical geometry is given by L z — L — NS, 
where 2S + 1 is the number of single-particle states on 
the sphere and 2S is chosen to be 2(A^ — 1) in order to 
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7/8 


3/4 


21/32 


19/32 


71/128 



TABLE I: The increase in energy of the excitation associated 
with the single-particle state <f> m (jn > 3) when the system 
develops from the nonvortex condensate with L — into the 
single- vortex condensate with L — N. 
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FIG. 2: Regular rotational bands above the single- vortex con- 
densate are shown, where the energy is measured from the line 
E = E L = N -5gN(L-N)/32 in units of gN . Among the rota- 
tional bands, three bands are shown by linking the collective 
excitations, Q\\L — TV) (solid curve), Q\\L — TV — 1) (dashed 
curve), and Q\\L — TV — 2) (dotted curve), for A < 7. The low- 
lying excitations can be understood as (multiple) excitations 
of these collective modes. Except for the states Qi\L — TV), 
Qi \L — TV — 1), and Q\ \L — TV — 2) (shown by squares), those 
states that include excitations of the center-of-mass motion 
are omitted. 




L 

FIG. 3: The excitation spectrum for < L < 56 with TV = 8, 
calculated in the spherical geometry. As the zero-interaction- 
energy state (i.e., the lowest energy state), the v — 1/2 Laugh- 
lin state appears at L — TV(TV — 1) = 56. In the inset, the ex- 
citation spectrum above the Laughlin state is shown against 
the total angular momentum in the spherical system. The 
lowest-energy branch exhibits a rotonlike structure without 
having a phonon branch despite the short-range nature of the 
interaction. 



study the v —1/2 Laughlin state. 

Figure [3] shows the excitation spectrum calculated in 
the spherical geometry for < L < 56 with N = 8. 
The v — 1/2 Laughlin state appears at L = TV(TV — 
1) = 56 (i.e., L z = 0) as the zero-interaction-energy state 
(i.e., the lowest energy state) for the contact repulsive 
interaction (Airg/R 2 ) J2i<j ^j)> where fli is a unit 

vector that specifies the location of the i— th boson on 
the sphere and the radius of the sphere is given by R = 
\J 5/2. The total AM, L to t > of this state is zero and the 
eigenstates with L to t = are plotted for L = NS. The 
eigenstates with non-zero total AM are plotted only for 
L = NS — Ltot without showing their (2Ltot + l)-fold 
degeneracy. 

The inset of Fig.|3shows the excitation spectrum above 
the Laughlin state as a function of L to t in the spheri- 
cal system. A rotonlike structure manifests itself as the 
lowest-energy excitation branch. We have confirmed that 
the excitation energy remains finite (1.229 ± 0.038 in 
units of g) in the thermodynamic limit |20| . We have 
also confirmed that the feature of the rotonlike minimum 
becomes more pronounced if we use long-range interac- 
tions such as the Coulomb interaction, Yli<j 
The two-roton continuum is also seen to exist above the 
single-roton branch. 

To understand the physics underlying these gap- 



ful excitations, let us recall that it is possible to 
produce quasi-particle excitations upon the Laughlin 
state 18] . One of them is the quasi- hole (vortex) ex- 
citation, rip[ e ~ |Zp|2/2 ( z p " z o)} Ili<j( z * - z j) 2 > and an- 
other one is the 'quasi-boson' (anti-vortex) excitation, 
ry e -M 2 /2 (J__ z *)] H^izi-Zj) 2 , where z denotes 
the location of the vortex or that of the anti- vortex. Then 
the lowest-lying neutral excitations are considered to be 
the bound states of a vortex-antivortex pair. 

By investigating a system with 25 = 2A^ — 3 or 
25 = 2iV— 1, we find that in the thermodynamic limit the 
quasi-boson creation energy remains finite (1.248 ±0.007 
in units of g) 20], while the quasi-hole creation energy 
is zero for repulsive contact interaction. Since each of 
quasiboson and quasihole excitations has AM N/2, the 
branch composed of a quasiboson-quasihole pair is ex- 
pected to end at L to t = N. This is indeed the case in 
the lowest-lying excitation branch as shown in Fig. 4, 
where the branches for system size N — 4, 5,6, 7 and 8 
are plotted against "wave number" k = Ltot/ R- 

In Fig. 4 it is seen that the minimum of the neutral- 
excitation energy (indicated by closed symbols) is real- 
ized at fc m ; n ~ 3 (Mui/h) 1 / 2 . For k > fe m i n , the excita- 
tion energy increases with increasing k and approaches 
the quasi-boson creation energy (indicated on the right 
side of the figure by a line for each N). As shown in the 
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FIG. 4: The lowest-lying rotonlike excitation branch is shown 
against "wave number" k = L to t/R for N = 4,5,6,7 and 
8. The minimum of the neutral-excitation energy is realized 
at fc m i n ~ 3 (Mui/h) 1 / 2 (indicated by closed symbols), and 
each branch ends at L to t = N as expected for a quasihole- 
quasiboson pair excitation. For k > k m i u , the excitation 
energy increases with increasing k and approaches the quasi- 
boson creation energy (indicated on the right side of the figure 
by a line for each N). In the inset, the quasiboson creation 
energy and the minimum value of the neutral-excitation en- 
ergy are shown against 1/(JV — 1) by open and closed circles, 
respectively. 

inset of Fig. 4, the difference between the quasi-boson 
creation energy (indicated by open circles) and the mini- 
mum value of the neutral-excitation energy (indicated by 
closed circles) decreases as the system size increases and 
may vanish in the thermodynamic limit. In the presence 
of this finite energy gap, however, this rotonlike minimum 
can survive even for very large values of N by introduc- 
ing, e.g., laser-induced dipole-dipole interactions between 
bosons |2lj . 

It should be noted that the lowest-lying spectrum has 
no phonon branch, despite the short-range nature of the 
interaction A similar roton spectrum has been re- 
ported in the fractional quantum Hall system |23| , where 
absence of the phonon branch has been attributed to the 
Higgs mechanism caused by the long-range nature of the 
Coulomb interaction. However, the rotonlike structure 
with no phonon branch found in this Letter for the con- 
tact interaction does not fit into this category. In the 
composite-boson picture, the v = 1/2 Laughlin state is 
nothing but a condensate of composite bosons, each of 
which is composed of a boson and two flux quanta at- 
tached to it. In the composite-particle picture for the 
quantum Hall system, it has been discussed that inter- 
actions between composite particles are different from 
those between original particles It might be con- 

ceived that the finite gap in the long-wavelength limit 



is related with momentum-dependent long-range interac- 
tions introduced on the composite-boson transformation 
(without being induced by a laser [HI)- However, the 
understanding of the gap formation found in this Letter 
remains to be clarified and merits further investigation. 
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